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Abstract 

Electron transfer is studied using a multi-level system coupled to a bosonic bath. Two- 
body correlation functions are obtained using both exact enumeration of spin paths and 
Monte Carlo simulation. We find that the phase boundary for the coherent-incoherent tran- 
sition lies at a smaller friction in the asymmetric two-level model than in the symmetric 
two-level model. A similar coherent-incoherent transition is observed for three-level system. 


1 Introduction 

Electron transfer in liquids is an important phenomenon in chemistry and physics. Following 
Marcus 5 picture [1], we use the spin boson model, which is a multi-level system coupled to a 
harmonic bath. The case of symmetric two-level model has been studied by many people, in the 
context of the Kondo problem and electron transfer in liquids[2]. It is not an exactly solvable 
model, except for the case of an adiabatic bath. In Section 3, we analyze the model with Feynman 
path integral which we evaluate with the exact enumeration of spin paths and Monte Carlo 
simulation. In Section 4, we present the results for an asymmetric two-level system and compare 
the results with the symmetric two-level system [3, 4, 5]. We also discuss the results for a three- 
level system. Finally, in Section 5, we summarize our results. 

2 Model Hamiltonian 

The Hamiltonian which defines our system is given by: 

H = H 0 + H b + H, nt . (1) 
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Here Ho is the the Hamiltonian for the free three-level system, 


Ho = 


1 Ei J\2 *Aa\ 
J\ 2 E 2 J23 
\ J \3 J 23 E 3 ) 


( 2 ) 


where Ei gives the energy for an electron localized on site i, and gives the electronic coupling 
between sites i and j. H B is the Hamiltonian of the harmonic bath. H tnt is the coupling 
between the three-level system and the bath, and is given by 


Hi n t = £12^12 + £23^23 ' t " £i 3 &i 3 - 
Here cr k ( k = (12), (23), (13)) is defined by 
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and the field £ k is a linear combination of the bath modes, £ k = JZjC^Xj. Harmonic baths linearly 
coupled to spin systems can be defined by the spectral density: 


c*c( 


l " rrijUj 


(5) 


Here, k and l are the pairs, (12), (23) or (13). For electron transfer in a liquid, an Ohmic spectral 
density gives a good model for the environment surrounding the electron. [2, 4] 


J h, (u ;) = T] kl ujexp{-ui/u^) 


( 6 ) 


3 Numerical Methods 

The quantity we are interested in is the following two-body time correlation function. 

< n!(0)n!(<) >= ]-Tr exp(— ftH)^ exp(i/ft/ft)n 1 exp(— iHt/h), 

z 


( 7 ) 


where nj is the population operator for being on site 1. Since H{ nt does not commute with Hq 
and H b , we use the following Suzuki- Trotter formula to do the path integral evaluations. 


exp (~/3H) = lim {exp{-(3H 0 /p) exp(-0H B /p)exp(-0H int /p)} p . 


p — *00 


( 8 ) 


Expanding this equation in a path integral representation and performing the Gaussian bath 
integrals analytically leads to: 


< ni(0)ni(t) >= ^ J2 ex P( ( f({ CT t}))ni(o-H-i) n i( a p+?+i)’ 


(9) 
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where p and q , sometimes called Trotter numbers, are the numbers of partitions the path has been 
divided into the thermal and real time part of the path. The sum goes over 3 A states of S=1 Ising 
system, where N is the total Trotter number, N = p+ 2 q. The action in Eq. 9 consists of two 
parts, 

V = 'r>Q + <fiint- 0 °) 

The first part, tp 0 , comes from the free three-level system Hamiltonian, and second part. q>, nt , 
results from the interaction between the system and the bath which we have already integrated 
out. The action ipi„ t is given by: 

Vm = (ll) 

1 ij kl 

where /*(cr*) is a quadratic function of a\ and the \ij is expressed by the spectral density J kl (u>) 
as follows. 

xS = - /“ J*V) cosh£(A, v - A,,). (12) 

7T JO l 

Here, 

A ,; = !>*, (13) 

fc=t 

where A* = fijp, —it/q or it/q, depending on the location of cr k on the path. 

Thus, the original three-level system coupled to a bath is now transformed into a S=1 Ising 
model with infinite- range interactions. This model is exactly solvable for the case of an adiabatic 
bath[6]. This case corresponds to Ising magnets in a very slow Gaussian field. In the general case, 
the analytic solution of this model is unknown. 

To do the numerical calculation on this model for small Trotter number, we used the exact 
enumeration of the path integral by adding all the 3^ states of the Ising spin system. For example, 
exact enumeration of A T =17 spins (1.29 x 10 8 states) takes 3 minutes on a Cray X-MP. 

We also used Monte Carlo simulation. Some filtering method[5] was necessary to overcome 
the so-called sign-problems [7, 8, 9, 10, 11, 12], which is often seen in quantum calculations of 
fermionic systems or spin models. 

4 Results 

To check the validity of our calculations, we have studied a free three-level system (i.e., coupling 
between the spin system and the bath is turned off), and a three-level system coupled to an 
adiabatic bath. For these systems, we computed the following correlation functions. 

Ci(f) = Re < ni(0)ni(f) >,Ci(t) = Re < ni(0)n 2 (f) >,Ca(f) = Re < ni(0)n 3 (f) > (14) 

Our results give good agreement with analytic results in these limits. 

For the special case of an asymmetric two-level system, we have calculated the time-correlation 
functions, and compared the results with the symmetric case[4]. Figure 1 shows the results 
for an asymmetric case with the parameters Ei = 0, E 2 = 2/i, E 3 = oc, and = — K , 
= ^23 = 0. We find from Figure 1 that the coherent-incoherent transition occurs around 
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r]/h = 0.4, or q = 2ri/fnr = 0.25. For a symmetric two-level system[4], the phase boundary for 
this temperature was at a = 0.4. Thus, the whole phase boundary is expected to lie at smaller 
q for the asymmetric case than for the symmetric case. In the asymmetric case, the symmetry is 
already broken, thus the coherence is easier to break than in the symmetric case. 



t 

FIG. 1 The correlation function C\(t) for an asymmetric two- level model of E 2 = 
2 1<, (3I< = 2.5, with rj/h =0.3 (circles), 0.4(squares), 0.5 (triangles) obtained by exact 
enumeration of spin paths with p=2 and q=7. 



FIG. 2 The correlation functions Ci(t) (circles), C 2 (t) (squares) and C 3 (t) (trian- 
gles) for a three-level system in the coherent region, J\ 2 = J 23 = —2K , J13 = —K, 

E x = 0.5 1<, E 2 = I<, E 3 = 0, PI< = 0.25, t l c 3 = 2, t] 2 = t 2 c 3 = 0, (t c = I</u> c ), 

T) 13 /h = i,Ti 12 /* = v 23 /b = 0- 

Figures 2 and 3 show results for a three-level system. The calculations were done by Monte 
Carlo simulation, with the Trotter numbers p=2 and q=10. As has been found in quantum Monte 
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Carlo simulation of spin systems[8], for systems having sign problems, it is usually more efficient 
to define 1 Monte Carlo step (MCS) to be a small subset of all the possible flips than to define it 
as all the possible flips. In our simulation, we define 1 MCS as 1 single spin-flip, 1 double spin- 
flip, 1 global spin-flip (i.e., flips all the spins), and 1 spin-flip of random length. We determine 
whether the spins should be flipped by the standard Metropolis algorithm, using the modulus 
of the complex weight exp (9?) for the transition probability. In this way, we have carried out 
simulation of 10 6 MCS, taking about 14 minutes on the Cray X-MP. To estimate the degree of 
sign cancellation, we measured the quantity r, the remaining ratio (related to the negative ratio 
defined in [8]). 


Z + -Z- 

T — 

z+ + z_ 


(15) 


Here Z+ denotes the sampled sum of the positive real parts of the weights, and Z_ denotes the 
same for absolute values of negative real parts of the weights. If r is small, the cancellation of the 
signs is large, leading to inaccuracy in the data. If r is large, the cancellation is small, thus giving 
more accurate results. In this definition, we are ignoring the effect of the cancellation due to the 
imaginary parts of the weights. 
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FIG. 3 The correlation functions Cj(t), Ci(t) and C 3 (f) for a three-level system in 
the incoherent region. J : , Ej and /? are same as in Fig. 2. t\ 3 = 2, t 12 = t 23 = 1, 

T} 13 /fl = 1, T] 12 /fl = r] 23 /h = 2. 

For the free three-level Hamiltonian H 0 , we assume the parameters, Ju = J23 = — 2 K , J 13 = — K, 
Ei = 0.5 K, E? = K, E 3 = 0. This could correspond to a system of redox-sites 1, 2 and 3, where 
distance between states 1 and 2, or 2 and 3 is shorter than the distance between 1 and 3. The 
correlation function < nj(0)n2(f) > approximately tells the rate of the electron transfer starting 
at state 1 and reaching 2 after time t. The energy of the state 2, E 2 , is assumed to be highest, 
followed by the energy of the state 1, E\. Starting from the state 1, the electron moves to state 2, 
since the exchange J 12 is stronger than J13, then gradually goes to state 3. In Fig. 3, the bath has 
the role of dephasing the coherence, preventing the electron population from going back to the 
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original state. This is a very brief picture of electron transfer over 3 states, with the intermediate 
state strongly coupled to the initial and terminal states. 

As for the effect of the sign cancellations, the remaining ratio r defined in Eq. 15 is 3% for 
Fig.2 and 13% for Fig.3. The magnitude of error is about 0.1 in Fig.2 and 0.02 in Fig. 3 The 
incoherent case has less effect of the exchange A , thus leading to less sign cancellations. 


5 Summary 

We have briefly described the numerical calculations of the time-correlation functions of an asym- 
metric two-level system and a three-level system. For an asymmetric two-level system, we find that 
the coherent-incoherent transition occurs at smaller friction 77 than for the symmetric case. For 
a three-level system, we calculated the population transfer of the elec tron when there is an inter- 
mediate high-energy state. We observed a coherent-incoherent transition similar to the two-level 
system. Further application of this model will be discussed elsewhere, 
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VI. Group Representations 
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